GWAS and meta-analysis identifies 49 genetic variants underlying critical COVID-19.

Critical illness in COVID-19 is an extreme and clinically homogeneous disease phenotype that we have previously shown1 to be highly efficient for discovery of genetic associations2. Despite the advanced stage of illness at presentation, we have shown that host genetics in patients who are critically ill with COVID-19 can identify immunomodulatory therapies with strong beneficial effects in this group3. Here we analyse 24,202 cases of COVID-19 with critical illness comprising a combination of microarray genotype and whole-genome sequencing data from cases of critical illness in the international GenOMICC (11,440 cases) study, combined with other studies recruiting hospitalized patients with a strong focus on severe and critical disease: ISARIC4C (676 cases) and the SCOURGE consortium (5,934 cases). To put these results in the context of existing work, we conduct a meta-analysis of the new GenOMICC genome-wide association study (GWAS) results with previously published data. We find 49 genome-wide significant associations, of which 16 have not been reported previously. To investigate the therapeutic implications of these findings, we infer the structural consequences of protein-coding variants, and combine our GWAS results with gene expression data using a monocyte transcriptome-wide association study (TWAS) model, as well as gene and protein expression using Mendelian randomization. We identify potentially druggable targets in multiple systems, including inflammatory signalling (JAK1), monocyte-macrophage activation and endothelial permeability (PDE4A), immunometabolism (SLC2A5 and AK5), and host factors required for viral entry and replication (TMPRSS2 and RAB2A).

Genetic predisposition to subjective well-being, depression, and suicide in relation to COVID-19 susceptibility and severity.

BACKGROUND: Epidemiological studies have reported associations between subjective well-being (SWB), depression, and suicide with COVID-19 illness, but the causality has not been established. We performed a two-sample Mendelian randomization (MR) analysis to investigate the causal link between SWB, depression, suicide and COVID-19 susceptibility and severity. METHODS: Summary statistics for SWB (298,420 cases), depression (113,769 cases) and suicide (52,208 cases) were obtained from three large-scale GWAS. Data on the associations between the Single Nucleotide Polymorphisms (SNPs) and COVID-19 (159,840 cases), hospitalized COVID-19 (44,986 cases), and severe COVID-19 (18,152 cases) were collected from the COVID-19 host genetics initiative. The causal estimate was calculated by the Inverse Variance Weighted, MR Egger and Weighted Median methods. Sensitivity tests were used to evaluate the validity of the causal relationship. RESULTS: Our results showed that genetically predicted SWB (OR = 0.98, 95 % CI: 0.86-1.10, P = 0.69), depression (OR = 0.76, 95 % CI: 0.54-1.06, P = 0.11), and suicide (OR = 0.99, 95 % CI: 0.96-1.02, P = 0.56) were not causally related to COVID-19 susceptibility. Similarly, we did not find a potential causal relationship between SWB, depression, suicide and COVID-19 severity. CONCLUSIONS: This indicated that positive or negative emotions would not make COVID-19 better or worse, and strategies that attempted to use positive emotions to improve COVID-19 symptoms may be useless. Improving knowledge about the SARS-CoV-2 and timely medical intervention to reduce panic during a pandemic is one of the effective measures to deal with the current decrease in well-being and increase in depression and suicide rates.

Associations between brain gene expression perturbations implicated by COVID-19 and psychiatric disorders.

BACKGROUND: Currently, there is increasing evidence from clinic, epidemiology, as well as neuroimaging, demonstrating neuropsychiatric abnormalities in COVID-19, however, whether there were associations between brain changes caused by COVID-19 and genetic susceptibility of psychiatric disorders was still unknown. METHODS: In this study, we performed a meta-analysis to investigate these associations by combing single-cell RNA sequencing datasets of brain tissues of COVID-19 and genome-wide association study summary statistics of psychiatric disorders. RESULTS: The analysis demonstrated that among ten psychiatric disorders, gene expression perturbations implicated by COVID-19 in excitatory neurons of choroid plexus were significantly associated with schizophrenia. CONCLUSIONS: Our analysis might provide insights for the underlying mechanism of the psychiatric consequence of COVID-19.

A MUC5B Gene Polymorphism, rs35705950-T, Confers Protective Effects Against COVID-19 Hospitalization but Not Severe Disease or Mortality.

Rationale: A common MUC5B gene polymorphism, rs35705950-T, is associated with idiopathic pulmonary fibrosis (IPF), but its role in severe acute respiratory syndrome coronavirus 2 infection and disease severity is unclear. Objectives: To assess whether rs35705950-T confers differential risk for clinical outcomes associated with coronavirus disease (COVID-19) infection among participants in the Million Veteran Program (MVP). Methods: The MUC5B rs35705950-T allele was directly genotyped among MVP participants; clinical events and comorbidities were extracted from the electronic health records. Associations between the incidence or severity of COVID-19 and rs35705950-T were analyzed within each ancestry group in the MVP followed by transancestry meta-analysis. Replication and joint meta-analysis were conducted using summary statistics from the COVID-19 Host Genetics Initiative (HGI). Sensitivity analyses with adjustment for additional covariates (body mass index, Charlson comorbidity index, smoking, asbestosis, rheumatoid arthritis with interstitial lung disease, and IPF) and associations with post-COVID-19 pneumonia were performed in MVP subjects. Measurements and Main Results: The rs35705950-T allele was associated with fewer COVID-19 hospitalizations in transancestry meta-analyses within the MVP (Ncases = 4,325; Ncontrols = 507,640; OR = 0.89 [0.82-0.97]; P = 6.86 × 10-3) and joint meta-analyses with the HGI (Ncases = 13,320; Ncontrols = 1,508,841; OR, 0.90 [0.86-0.95]; P = 8.99 × 10-5). The rs35705950-T allele was not associated with reduced COVID-19 positivity in transancestry meta-analysis within the MVP (Ncases = 19,168/Ncontrols = 492,854; OR, 0.98 [0.95-1.01]; P = 0.06) but was nominally significant (P < 0.05) in the joint meta-analysis with the HGI (Ncases = 44,820; Ncontrols = 1,775,827; OR, 0.97 [0.95-1.00]; P = 0.03). Associations were not observed with severe outcomes or mortality. Among individuals of European ancestry in the MVP, rs35705950-T was associated with fewer post-COVID-19 pneumonia events (OR, 0.82 [0.72-0.93]; P = 0.001). Conclusions: The MUC5B variant rs35705950-T may confer protection in COVID-19 hospitalizations.

SUMMIT: An integrative approach for better transcriptomic data imputation improves causal gene identification.

Genes with moderate to low expression heritability may explain a large proportion of complex trait etiology, but such genes cannot be sufficiently captured in conventional transcriptome-wide association studies (TWASs), partly due to the relatively small available reference datasets for developing expression genetic prediction models to capture the moderate to low genetically regulated components of gene expression. Here, we introduce a method, the Summary-level Unified Method for Modeling Integrated Transcriptome (SUMMIT), to improve the expression prediction model accuracy and the power of TWAS by using a large expression quantitative trait loci (eQTL) summary-level dataset. We apply SUMMIT to the eQTL summary-level data provided by the eQTLGen consortium. Through simulation studies and analyses of genome-wide association study summary statistics for 24 complex traits, we show that SUMMIT improves the accuracy of expression prediction in blood, successfully builds expression prediction models for genes with low expression heritability, and achieves higher statistical power than several benchmark methods. Finally, we conduct a case study of COVID-19 severity with SUMMIT and identify 11 likely causal genes associated with COVID-19 severity.

Inflammasome Genetic Variants Are Associated with Protection to Clinical Severity of COVID-19 among Patients from Rio de Janeiro, Brazil.

COVID-19 has a broad spectrum of clinical manifestations, from asymptomatic or mild/moderate symptoms to severe symptoms and death. The mechanisms underlying its clinical evolution are still unclear. Upon SARS-CoV-2 infection, host factors, such as the inflammasome system, are activated by the presence of the virus inside host cells. The search for COVID-19 risk factors is of relevance for clinical management. In this study, we investigated the impact of inflammasome single-nucleotide polymorphisms (SNPs) in SARS-CoV-2-infected individuals with distinct severity profiles at clinical presentation. Patients were divided into two groups according to disease severity at clinical presentation based on the WHO Clinical Progression Scale. Group 1 included patients with mild/moderate disease (WHO < 6; n = 76), and group 2 included patients with severe/critical COVID-19 (WHO ≥ 6; n = 357). Inpatients with moderate to severe/critical profiles were recruited and followed-up at Hospital Center for COVID-19 Pandemic - National Institute of Infectology (INI)/FIOCRUZ, RJ, Brazil, from June 2020 to March 2021. Patients with mild disease were recruited at Oswaldo Cruz Institute (IOC)/FIOCRUZ, RJ, Brazil, in August 2020. Genotyping of 11 inflammasome SNPs was determined by real-time PCR. Protection and risk estimation were performed using unconditional logistic regression models. Significant differences in NLRP3 rs1539019 and CARD8 rs2043211 were observed between the two groups. Protection against disease severity was associated with the A/A genotype (ORadj = 0.36; P = 0.032), allele A (ORadj = 0.93; P = 0.010), or carrier-A (ORadj = 0.45; P = 0.027) in the NLRP3 rs1539019 polymorphism; A/T genotype (ORadj = 0.5; P = 0.045), allele T (ORadj = 0.93; P = 0.018), or carrier-T (ORadj = 0.48; P = 0.029) in the CARD8 rs2043211 polymorphism; and the A-C-G-C-C (ORadj = 0.11; P = 0.018), A-C-G-C-G (ORadj = 0.23; P = 0.003), C-C-G-C-C (ORadj = 0.37; P = 0.021), and C-T-G-A-C (ORadj = 0.04; P = 0.0473) in NLRP3 genetic haplotype variants. No significant associations were observed for the other polymorphisms. To the best of our knowledge, this is the first study demonstrating an association between CARD8 and NLRP3 inflammasome genetic variants and protection against COVID-19 severity, contributing to the discussion of the impact of inflammasomes on COVID-19 outcomes.

A genetic variant in IL-6 lowering its expression is protective for critical patients with COVID-19.

Critical coronavirus disease 2019 (COVID-19) is associated with high mortality and potential genetic factors have been reported to be involved in the development of critical COVID-19. We performed a genome-wide association study to identify the genetic factors responsible for developing critical COVID-19. 632 critical patients with COVID-19 and 3021 healthy controls from the Chinese population were recruited. First, we identified a genome-wide significant difference of IL-6 rs2069837 (p = 9.73 × 10-15, OR = 0.41) between 437 critical patients with COVID-19 and 2551 normal controls in the discovery cohort. When replicated these findings in a set of 195 patients with critical COVID-19 and 470 healthy controls, we detected significant association of rs2069837 with COVID-19 (p = 8.89 × 10-3, OR = 0.67). This variant surpassed the formal threshold for genome-wide significance (combined p = 4.64 × 10-16, OR = 0.49). Further analysis revealed that there was a significantly stronger expression of IL-6 in the serum from patients with critical COVID-19 than in that from patients with asymptomatic COVID-19. An in vitro assay showed that the A to G allele changes in rs2069837 within IL-6 obviously decreased the luciferase expression activity. When analyzing the effect of this variant on the IL-6 in the serum based on the rs2069837 genotype, we found that the A to G variation in rs2069837 decreased the expression of IL-6, especially in the male. Overall, we identified a genetic variant in IL-6 that protects against critical conditions with COVID-19 though decreasing IL-6 expression in the serum.

Association of polymorphisms in tumor necrosis factors with SARS-CoV-2 infection and mortality rate: A case-control study and in silico analyses.

The present coronavirus disease 2019 (COVID-19) is spreading rapidly and existing data has suggested a number of susceptibility factors for developing a severe course of the disease.  The current case-control experiment is aimed to study the associations of genetic polymorphisms in tumor necrosis factors (TNFs) with COVID-19 and its mortality rate. A total of 550 participants (275 subjects and 275 controls) were enrolled. The tetra-amplification refractory mutation system polymerase chain reaction technique was recruited to detect -308G>A TNFα and +252A>G TNFß polymorphisms among the Iranian subjects. We demonstrated that carriers of the G allele of TNFß-252A/G, rs909253 A>G were more frequent in COVID-19 subjects compared to the healthy group and this allele statistically increased the disease risk (odds ratio [OR] = 1.55, 95% confidence interval [CI] = 1.23-1.96, p < 0.0001). At the same time, the A allele of TNFα-311A/G, rs1800629 G>A moderately decreased the risk of COVID-19 (OR = 0.68, 95% CI = 0.53-0.86, p < 0.002). Also, we analyzed the various genotypes regarding the para-clinical and disorder severity; we found that in the AA genotype of TNFß-252A/G (rs909253 A>G), the computed tomography scan pattern was different in comparison to cases carrying the AG genotype with p1 < 0.001. In addition, in the severe cases of COVID-19, leukocyte and neutrophil count and duration of intensive care unit hospitalization in the deceased patients were significantly increased (p < 0.001). Moreover, the TNFα-311A/G (rs1800629 G>A) variant is likely to change the pattern of splicing factor sites. Our findings provided deep insights into the relationship between TNFα/TNFß polymorphisms and severe acute respiratory syndrome coronavirus 2. Replicated studies may give scientific evidence for exploring molecular mechanisms of COVID-19 in other ethnicities.

Genetic Predisposition and the Variable Course of Infectious Diseases.

BACKGROUND: Contact with a pathogen is followed by variable courses of infectious disease, which are only partly explicable by classical risk factors. The susceptibility to infection is variable, as is the course of disease after infection. In this review, we discuss the extent to which this variation is due to genetic factors of the affected individual (the host). METHODS: Selective review of the literature on host genetics in infectious disease, with special attention to the pathogens SARSCoV- 2, influenza viruses, Mycobacterium tuberculosis, and human immunodeficiency virus (HIV). RESULTS: Genetic variants of the host contribute to the pathogenesis of infectious diseases. For example, in HIV infection, a relatively common variant leading to a loss of function of the HIV co-receptor CCR5 affects the course of the disease, as do variants in genes of the major histocompatibility complex (MHC) region. Rare monogenic variants of the interferon immune response system contribute to severe disease courses in COVID-19 and influenza (type I interferon in these two cases) and in tuberculosis (type II interferon). An estimated 1.8% of life-threatening courses of COVID-19 in men under age 60 are caused by a deficiency of toll-like receptor 7. The scientific understanding of host genetic factors has already been beneficial to the development of effective drugs. In a small number of cases, genetic information has also been used for individual therapeutic decision-making and for the identification of persons at elevated risk. CONCLUSION: A comprehensive understanding of host genetics can improve the care of patients with infectious diseases. Until the present, the clinical utility of host genetics has been limited to rare cases; in the future, polygenic risk scores summarizing the relevant genetic variants in each patient will enable a wider benefit. To make this possible, multicenter studies are needed that will systematically integrate clinical and genetic data.

SARS-CoV-2-Encoded MiRNAs Inhibit Host Type I Interferon Pathway and Mediate Allelic Differential Expression of Susceptible Gene.

Infection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing the rapid spread of coronavirus disease 2019 (COVID-19), has generated a public health crisis worldwide. The molecular mechanisms of SARS-CoV-2 infection and virus-host interactions are still unclear. In this study, we identified four unique microRNA-like small RNAs encoded by SARS-CoV-2. SCV2-miR-ORF1ab-1-3p and SCV2-miR-ORF1ab-2-5p play an important role in evasion of type I interferon response through targeting several genes in type I interferon signaling pathway. Particularly worth mentioning is that highly expressed SCV2-miR-ORF1ab-2-5p inhibits some key genes in the host innate immune response, such as IRF7, IRF9, STAT2, OAS1, and OAS2. SCV2-miR-ORF1ab-2-5p has also been found to mediate allelic differential expression of COVID-19-susceptible gene OAS1. In conclusion, these results suggest that SARS-CoV-2 uses its miRNAs to evade the type I interferon response and links the functional viral sequence to the susceptible genetic background of the host.

Prognostic impact of toll-like receptors gene polymorphism on outcome of COVID-19 pneumonia: A case-control study.

Toll-like receptor 3 (TLR3) and TLR7 genes are involved in the host immune response against viral infections including SARS-COV-2. This study aimed to investigate the association between the TLR3(rs3775290) and TLR7(rs179008) polymorphisms with the prognosis and susceptibility to COVID-19 pneumonia accompanying SARS-COV-2 infection. This case-control study included 236 individuals: 136 COVID-19 pneumonia patients and 100 age and sex-matched controls. Two polymorphisms (TLR3 rs3775290 and TLR7 rs179008) were genotyped by allelic discrimination through TaqMan real-time PCR. This study also investigated predictors of mortality in COVID-19 pneumonia through logistic regression. The mutant 'T/T' genotypes and the 'T' alleles of TLR3(rs3775290) and TLR7(rs179008) polymorphisms were significantly associated with increased risk of COVID-19 pneumonia. This study did not report association between the mutant 'T/T' genotypes of TLR3(rs3775290) and TLR7(rs179008) and the disease outcome. In multivariate analysis, the independent predictors of mortality in COVID-19 pneumonia were male sex, SPO2 ≤ 82%, INR > 1, LDH ≥ 1000 U/l, and lymphocyte count<900/mm3 (P < 0.05).

Implications of the Immune Polymorphisms of the Host and the Genetic Variability of SARS-CoV-2 in the Development of COVID-19.

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is responsible for the current pandemic affecting almost all countries in the world. SARS-CoV-2 is the agent responsible for coronavirus disease 19 (COVID-19), which has claimed millions of lives around the world. In most patients, SARS-CoV-2 infection does not cause clinical signs. However, some infected people develop symptoms, which include loss of smell or taste, fever, dry cough, headache, severe pneumonia, as well as coagulation disorders. The aim of this work is to report genetic factors of SARS-CoV-2 and host-associated to severe COVID-19, placing special emphasis on the viral entry and molecules of the immune system involved with viral infection. Besides this, we analyze SARS-CoV-2 variants and their structural characteristics related to the binding to polymorphic angiotensin-converting enzyme type 2 (ACE2). Additionally, we also review other polymorphisms as well as some epigenetic factors involved in the immunopathogenesis of COVID-19. These factors and viral variability could explain the increment of infection rate and/or in the development of severe COVID-19.

Association of genetic variations in ACE2, TIRAP and factor X with outcomes in COVID-19.

BACKGROUND: Coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), can manifest with varying disease severity and mortality. Genetic predisposition influences the clinical course of infectious diseases. We investigated whether genetic polymorphisms in candidate genes ACE2, TIRAP, and factor X are associated with clinical outcomes in COVID-19. METHODS: We conducted a single-centre retrospective cohort study. All patients who visited the emergency department with SARS-CoV-2 infection proven by polymerase chain reaction were included. Single nucleotide polymorphisms in ACE2 (rs2285666), TIRAP (rs8177374) and factor X (rs3211783) were assessed. The outcomes were mortality, respiratory failure and venous thromboembolism. Respiratory failure was defined as the necessity of >5 litres/minute oxygen, high flow nasal oxygen suppletion or mechanical ventilation. RESULTS: Between March and April 2020, 116 patients (35% female, median age 65 [inter quartile range 55-75] years) were included and treated according to the then applicable guidelines. Sixteen patients (14%) died, 44 patients (38%) had respiratory failure of whom 23 required endotracheal intubation for mechanical ventilation, and 20 patients (17%) developed venous thromboembolism. The percentage of TIRAP polymorphism carriers in the survivor group was 28% as compared to 0% in the non-survivor group (p = 0.01, Bonferroni corrected p = 0.02). Genotype distribution of ACE2 and factor X did not differ between survivors and non-survivors. CONCLUSION: This study shows that carriage of TIRAP polymorphism rs8177374 could be associated with a significantly lower mortality in COVID-19. This TIRAP polymorphism may be an important predictor in the outcome of COVID-19.

Genome-wide association study of hospitalized COVID-19 patients in the United Arab Emirates.

BACKGROUND: The heterogeneity in symptomatology and phenotypic profile attributable to COVID-19 is widely unknown. The objective of this manuscript is to conduct a trans-ancestry genome wide association study (GWAS) meta-analysis of COVID-19 severity to improve the understanding of potentially causal targets for SARS-CoV-2. METHODS: This cross-sectional study recruited 646 participants in the UAE that were divided into two phenotypic groups based on the severity of COVID-19 phenotypes, hospitalized (n=482) and non-hospitalized (n=164) participants. Hospitalized participants were COVID-19 patients that developed acute respiratory distress syndrome (ARDS), pneumonia or progression to respiratory failure that required supplemental oxygen therapy or mechanical ventilation support or had severe complications such as septic shock or multi-organ failure. We conducted a trans-ancestry meta-analysis GWAS of European (n=302), American (n=102), South Asian (n=99), and East Asian (n=107) ancestry populations. We also carried out comprehensive post-GWAS analysis, including enrichment of SNP associations in tissues and cell-types, expression quantitative trait loci and differential expression analysis. FINDINGS: Eight genes demonstrated a strong association signal: VWA8 gene in locus 13p14·11 (SNP rs10507497; p=9·54 x10-7), PDE8B gene in locus 5q13·3 (SNP rs7715119; p=2·19 x10-6), CTSC gene in locus 11q14·2 (rs72953026; p=2·38 x10-6), THSD7B gene in locus 2q22·1 (rs7605851; p=3·07x10-6), STK39 gene in locus 2q24·3 (rs7595310; p=4·55 x10-6), FBXO34 gene in locus 14q22·3 (rs10140801; p=8·26 x10-6), RPL6P27 gene in locus 18p11·31 (rs11659676; p=8·88 x10-6), and METTL21C gene in locus 13q33·1 (rs599976; p=8·95 x10-6). The genes are expressed in the lung, associated to tumour progression, emphysema, airway obstruction, and surface tension within the lung, as well as an association to T-cell-mediated inflammation and the production of inflammatory cytokines. INTERPRETATION: We have discovered eight highly plausible genetic association with hospitalized cases in COVID-19. Further studies must be conducted on worldwide population genetics to facilitate the development of population specific therapeutics to mitigate this worldwide challenge. FUNDING: This review was commissioned as part of a project to study the host cell receptors of coronaviruses funded by Khalifa University's CPRA grant (Reference number 2020-004).

Association of KIR gene polymorphisms with COVID-19 disease.

BACKGROUND: Natural killer (NK) cells play an essential role against viruses. NK cells express killer cell immunoglobulin-like receptors (KIRs) which regulate their activity and function. The polymorphisms in KIR haplotypes confer differential viral susceptibility and disease severity caused by infections. We investigated the association between KIR genes and COVID-19 disease severity. METHODS: 424 COVID-19 positive patients were divided according to their disease severity into mild, moderate and severe. KIR genes were genotyped using next generation sequencing (NGS). Association between KIR genes and COVID-19 disease severity was conducted and significant correlations were reported. RESULTS: In the COVID-19 patients, KIR Bx genotype was more common than AA genotype. The Bx genotype was found more frequently in patients with mild disease, while in severe disease the AA genotype was more common than the Bx genotype. The KIR2DS4 gene carried the highest risk for severe COVID-19 infection (OR 8.48, pc= 0.0084) followed by KIR3DL1 (OR 7.61, pc= 0.0192). CONCLUSIONS: Our findings suggest that KIR2DS4 and KIR3DL1 genes carry risk for severe COVID-19 disease.

APOL1 risk variants in individuals of African genetic ancestry drive endothelial cell defects that exacerbate sepsis.

The incidence and severity of sepsis is higher among individuals of African versus European ancestry. We found that genetic risk variants (RVs) in the trypanolytic factor apolipoprotein L1 (APOL1), present only in individuals of African ancestry, were associated with increased sepsis incidence and severity. Serum APOL1 levels correlated with sepsis and COVID-19 severity, and single-cell sequencing in human kidneys revealed high expression of APOL1 in endothelial cells. Analysis of mice with endothelial-specific expression of RV APOL1 and in vitro studies demonstrated that RV APOL1 interfered with mitophagy, leading to cytosolic release of mitochondrial DNA and activation of the inflammasome (NLRP3) and the cytosolic nucleotide sensing pathways (STING). Genetic deletion or pharmacological inhibition of NLRP3 and STING protected mice from RV APOL1-induced permeability defects and proinflammatory endothelial changes in sepsis. Our studies identify the inflammasome and STING pathways as potential targets to reduce APOL1-associated health disparities in sepsis and COVID-19.

Poking COVID-19: Insights on Genomic Constraints among Immune-Related Genes between Qatari and Italian Populations.

Host genomic information, specifically genomic variations, may characterize susceptibility to disease and identify people with a higher risk of harm, leading to better targeting of care and vaccination. Italy was the epicentre for the spread of COVID-19 in Europe, the first country to go into a national lockdown and has one of the highest COVID-19 associated mortality rates. Qatar, on the other hand has a very low mortality rate. In this study, we compared whole-genome sequencing data of 14398 adults and Qatari-national to 925 Italian individuals. We also included in the comparison whole-exome sequence data from 189 Italian laboratory-confirmed COVID-19 cases. We focused our study on a curated list of 3619 candidate genes involved in innate immunity and host-pathogen interaction. Two population-gene metric scores, the Delta Singleton-Cohort variant score (DSC) and Sum Singleton-Cohort variant score (SSC), were applied to estimate the presence of selective constraints in the Qatari population and in the Italian cohorts. Results based on DSC and SSC metrics demonstrated a different selective pressure on three genes (MUC5AC, ABCA7, FLNA) between Qatari and Italian populations. This study highlighted the genetic differences between Qatari and Italian populations and identified a subset of genes involved in innate immunity and host-pathogen interaction.

COVID-19 genetic risk variants are associated with expression of multiple genes in diverse immune cell types.

Common genetic polymorphisms associated with COVID-19 illness can be utilized for discovering molecular pathways and cell types driving disease pathogenesis. Given the importance of immune cells in the pathogenesis of COVID-19 illness, here we assessed the effects of COVID-19-risk variants on gene expression in a wide range of immune cell types. Transcriptome-wide association study and colocalization analysis revealed putative causal genes and the specific immune cell types where gene expression is most influenced by COVID-19-risk variants. Notable examples include OAS1 in non-classical monocytes, DTX1 in B cells, IL10RB in NK cells, CXCR6 in follicular helper T cells, CCR9 in regulatory T cells and ARL17A in TH2 cells. By analysis of transposase accessible chromatin and H3K27ac-based chromatin-interaction maps of immune cell types, we prioritized potentially functional COVID-19-risk variants. Our study highlights the potential of COVID-19 genetic risk variants to impact the function of diverse immune cell types and influence severe disease manifestations.

SARS-CoV-2 susceptibility and COVID-19 disease severity are associated with genetic variants affecting gene expression in a variety of tissues.

Variability in SARS-CoV-2 susceptibility and COVID-19 disease severity between individuals is partly due to genetic factors. Here, we identify 4 genomic loci with suggestive associations for SARS-CoV-2 susceptibility and 19 for COVID-19 disease severity. Four of these 23 loci likely have an ethnicity-specific component. Genome-wide association study (GWAS) signals in 11 loci colocalize with expression quantitative trait loci (eQTLs) associated with the expression of 20 genes in 62 tissues/cell types (range: 1:43 tissues/gene), including lung, brain, heart, muscle, and skin as well as the digestive system and immune system. We perform genetic fine mapping to compute 99% credible SNP sets, which identify 10 GWAS loci that have eight or fewer SNPs in the credible set, including three loci with one single likely causal SNP. Our study suggests that the diverse symptoms and disease severity of COVID-19 observed between individuals is associated with variants across the genome, affecting gene expression levels in a wide variety of tissue types.